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IONOSPHERE TOPSIDE SOUNDER STUDIES, I: 
THE REDUCTION OF ALOUETTE I IONOGRAMS 


TO ELECTRON DENSITY PROFILES 

By J. O. Thomas, B. R. Briggs, L. Colin, 
M. J. Ryeroft, and Margaret Covert 
Ames Research Center 


SUMMARY 




In radar soundings of the topside of the ionosphere, an ionosonde in a 
vehicle immersed in the plasma above the peak of the F2 layer measures, as a 
function of frequency, f, the virtual depth, h' , below the vehicle, at which 
radio waves are reflected from the topside of the ionosphere. In this report 
a number of mathematical procedures and digital computer programs for convert 
ing the observed topside h'(f) curves into electron density profiles are 
described. 


INTRODUCTION 




The experiments carried on board the Canadian satellite Alouette have 
been described in reference 1 and by Warren (refs. 2 and 3)* The first 
results have been presented in a series of papers in the Canadian Journal of 
Physics (refs. 4 to 9)» Details concerning the construction, operation and 
orbit characteristics of the satellite can be found in the references quoted 
above . 

In the topside sounder experiment (see section entitled "The Basic Prob- 
lem") the virtual depths of reflection below the satellite, h' , of radio 
pulses vertically incident on the topside ionosphere are measured as a func- 
tion of frequency, f. The resulting topside h’(f) data, when displayed on 
photographic paper or film, are known as ionograms. These ionograms contain 
information about the distribution of electron density with height in the 
lower exosphere, a detailed knowledge of which is Important for the physics 
of the ionosphere. The derivation of the electron density distribution as a 
function of height (N(h) profiles) from the ionogram is, however, a complicated 
procedure best carried out by a digital computer program. 

This report presents information concerning such calculations which have 
been programmed at the NASA Ames Research Center for an IBM 7094 digital com- 
puter. These programs, in conjunction with a supplementary program described 
in a companion paper (ref. 10), are particularly suitable for the routine 
reduction of Alouette I ionograms to topside N(h) profiles and form an 
extension of work described recently by Thomas and Sader (refs. 11 and 12), 
Thomas, Long and Westover (ref. 13) and Thomas and Westover (ref. 14) . 


The basic problem is first described. The relevant calculations and J 
procedures used for testing the computations together with details concerning 
the operational procedures to be employed for converting an ionogram to an 
electron density profile are subsequently presented. 

In the last section the advantages and disadvantages of each technique 
are discussed, and test and sample cases presented in graphic and tabular 
form. The tabulated h ! (f) curves and corresponding N(h) profiles computed 
both for model and real ionospheres can be used for comparison purposes and 
for testing newly developed programs. 

Appendixes are included which prescribe standardized ionogram scaling or 
digitizing procedures and describe in detail the existing h r (f) input and 
N(h) output formats. 

The authors wish to acknowledge gratefully the courtesy of scientists of 
the Canadian Defence Research Telecommunications Establishment, Ottawa, Canada, 
particularly Dr. J. H. Chapman and Dr. G. L. Nelms who made a number of top- 
side ionograms available to the authors in the early stages of development of 
this work. They wish also to thank the staff of the Radioscience Laboratory, 
Stanford University, for recording ionograms. 


NOMENCLATURE 


The symbols listed below define physical quantities. Symbols introduced 
in mathematical derivations are defined as needed in the text and are not 
listed here. 

c speed of light 

f frequency of a radio pulse 

f 1 defined as - f x f H for f x > % 

f electron gyrofrequency 
fjj plasma frequency 

F total magnitude of the geomagnetic field 

G 0 defined as ^ Jl - (%/fi ) 2 

G x defined as - (%/q ) 2 

h real height above the surface of the earth 
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h' % virtual depth below the satellite 

H' group retardation 

I magnetic dip angle 

N electron density 

t time 

u group velocity 

Ah real depth measured positively downward from the vehicle 

Pq group refractive index of the Ordinary ray 

p^. group refractive index of the Extraordinary ray 

Subscripts 

i discrete value 

o quantity related to the Ordinary trace 

x quantity related to the Extraordinary trace 

v quantity measured at the position of the vehicle 

THE BASIC PROBLEM 


In ionospheric variable frequency vertical incidence radio-pulse sound- 
ings from the ground or from satellite -borne radars, the time delay t 
between the emission of a radar pulse and the reception of an echo after its 
round trip from the transmitter to the point of reflection and back is meas- 
ured. The virtual height of reflection above the ground, or, for a topside 
sounder, the virtual depth of reflection below the satellite h', is defined 

t>y 


h' = | ct (1) 

where c is the velocity of light. The virtual height or depth thus corres- 
ponds to the distance which the radio pulse would have covered if it had 
everywhere along its path traveled with the velocity of light. 

In a plasma, the real speed of the pulse is given by the group velocity 
u so that 
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where h is the real distance traveled by the pulse. Thus the total time of 
flight of the pulse and echo is given by 


so that 


ph dh 

J o u 


h 1 = c 


ph dh 

J 0 u 


(3 


(4 


The group refractive index p' is defined as c/u and can be calculated for 
the Ordinary and Extraordinary rays from the Appleton-Hartree magneto-ionic 
equation. 

Hence, showing the functional dependencies of the quantities concerned 
we may write, neglecting the effect of collisions and measuring h and h' 
from the ground. 


h'(f) = J 


h(N) 

O 


p'(f,N,I,F)dh 


(5 


where 

f frequency of radio pulse 
N electron density (a function of height) 

I magnetic dip angle 

F total magnitude of the earth 1 s magnetic field 

Equation (5) may also be related to the plasma frequency, f^, and elec- 
tron gyrofrequency, fg, by the relations 

1.2388x10 4 f^ and f H — 2.7994 F 

with fjj and fjj given in Mc/s, N in electrons per cc and F in gauss. 

Thus, 

h’(f) = J h(f N ) p-(f,f N ,I,f H )dh (6 



This is the basic equation relating the observed h’(f) data to the N(h) 
profile. In practice it is necessary to invert equation (6) so that h(fjj) 
can be determined given h T (f). In order to do this, it is convenient to 
transform (6) so that for the measurements of the Ordinary ray 

h°(fi) = (§~) df N (>) 

The subscript i indicates a value of frequency corresponding to a par- 
ticular value , h±, of real height and to a particular value, h^, of the vir- 
tual height measured on an ionogram. The upper limit of integration in 
equation (7a) is the plasma frequency in the region where the Ordinary ray of 
frequency fj_ Is reflected. For the Extraordinary ray 

t 

- J* 1 4 (fjj) 4% (To) 

with 

*1 = f ii - ( f xi > f Hi) (8) 


the subscript x being included to indicate that the frequency is now meas- 
ured on the Extraordinary trace. 


A particular frequency of fxi> equal to f x t, say measured on the Extraor- 
dinary trace, is reflected from the same real height as an Ordinary ray of 
frequency f a , say, if, and only if, 


.1 2 


The frequencies f a and f ^ 


= f 2 - f _.f. 


xb 


xb Hb 


_ f 2 _ f 2 

“ a N 


are then known as "corresponding" frequencies. 


In general, one can write (see eq. (6)) 

' 9 (%) 


■ I 


M-’ (f 


dh 

dcp 


dcp 


( 9 ) 


where cp is some single valued function of fjj, and fg is the gyrofrequency 
for electrons. Equations (7a) and (7b) obviously correspond to the case 


It should be noted that dh/df N must remain finite over the interval 
concerned so that for the inversion of the h'(f) curve it is required that 
h(fjr) or, generally, h(cp) be monotonic. 

Equations (7a) and (7b) apply for ground-based soundings. They may be 
rewritten, with slight notational changes, for application to topside iono- 
grams , as follows : 
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( 10 ) 



( 11 ) 


The subscript v in the lower limits of these two integrals indicates that 
the plasma frequency is that at the position of the vehicle, and Ah is used 
here for the real depth of reflection below the satellite. Equation ( 9 ) may 
be altered similarly giving, for the Ordinary ray, for example 


cp( f^) 

p;(f,f N ,i,f H ) f^Vcp (12) 

cp(f v ) V d V 


From equation (10), (ll) or (12) it is possible to deduce Ah(fjj) from the 
observed h’(f) variation recorded on the topside ionogram. This may sub- 
sequently be converted into an electron density profile giving N as a func- 
tion of h, the height above the ground by means of the relation Ah = h v - h 
where h v is the height of the satellite. 

In procedures described in the 
subsequent pages for the inversion of 
equations (10), (ll) and (12) it is 
assumed that the plasma frequency f v 
or f y at the vehicle is known, as 
are also the quantities I, fjj and h v . 
Techniques for determining these four 
parameters are described in refer - 
erence 10. It should be noted that in 
computing fj_ or f-J. from equation (8) 
the appropriate value fjjq or fjfv of 
the height -dependent gyrofrequency 
should be used. 

Figure 1 illustrates schematically 
a typical topside profile with associ- 
ated Ordinary and Extraordinary ray 
traces . 



Fig. 1. - Schematic diagram of a typical topside 
electron density profile [f n (Ah)] and its 
associated Extraordinary [hi(f)] and Ordi- 
nary [ho(f)] ionogram traces. In this 
diagram f xi has teen chosen as the "corre- 
sponding frequency to f-j. (see text). 
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DIGITAL COMPUTER PROGRAMS FOR THE REDUCTION OF 
TOPSIDE IONOGRAMS TO ELECTRON DENSITY PROFILES 


General Remarks 

Several programs are available at the NASA Ames Research Center for 
reducing topside ionogram data. The techniques used are noted below and are 
directly based on ; or are simple developments of, analytical procedures 
already extensively described in the scientific literature. The appropriate 
relevant references are given in table I for each technique. 

1. Single polynomial method 

2. Overlapping polynomial method 

3. Lamination method 

4. Contiguous parabola method 

The input data include the dip angle at the vehicle the gyrofrequency at 

the position of the satellite, fg v , and a number of pairs of frequency and 
virtual -depth data as read from an ionogram according to the reading or scal- 
ing procedure described in appendix A. Identification and control information 
is also read in. The input data are punched onto cards, in a single format 
which can be used for all the programs currently available. The punched card 
input formats are described in detail In appendix E. 

The basic output for a single ionogram includes all of the input quan- 
tities, and the calculated real depth and electron density. Options are avail- 
able in the overlapping polynomial program to print interpolated summary 
tables of electron density and scale height, H, as functions of real depth at 
50-km intervals for a series of ionograms. The scale height H is approxi- 
mated by |N/(AN/Ah)|. Details of output formats for the several programs, as 
well as computer time estimates, are given In appendix 0. 


The Single Polynomial Method 

The topside single polynomial program is essentially identical to that 
described briefly by Thomas, Long and Westover (ref.13) and in more detail in 
a report by Thomas and Westover (ref. 14). Basically the program consists of 
a simple development of the single polynomial technique (for references see 
table I) proposed by TItheridge (ref. 15) for the reduction of bottoms ide 
ionograms. The modification outlined briefly below allows for the two con- 
ditions appropriate to the topside case, namely, that both N and dN/dh are 
finite at Ah = 0. 

It is assumed that the N(h) curve (or part of the N(h)) curve can be 
approximated by the polynomial 
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Ah = 


n 


l 


a^f. 


N 


-Nv' 


j 

(13) 


where fjj v is the plasma frequency at the vehicle and n is the number of 
points through which the polynomial is to be fitted. Details are omitted here, 
but it is shown in the references noted above that Ah can be determined by 
solution of the matrix equation 


[Ah] = [A][B _1 ][h«] (14) 


The elements of [A] are 

a ij = ( f i " f Nv) * (15) 

For the Ordinary ray the elements of [ B] are 

r> f i 

= j / - (l6) 

J f Nv 


Equation (l6) may be written in the form 

pit / 2 


^ij - ^i 


sin" 1 (f Nv /f i ) 


G 0 (f“i sin 6 - f Nv ) 


J-i 


d0 


(17) 


where the transformation fjj = fj_ sin 0 has been made, and where Pq has 
been replaced by the function G 0 /VT ' ( f |/ f i) • 

The above equations apply for the Extraordinary trace if fj_, p^ and G Q 
are replaced by p^ and G x , respectively, where 


and 


f 1 2 - f 2 . - f . f 
I i " x xi I xi I Hv 


Gi 


P-x- = 


Jl - (fg/fi 2 ) 
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|j.q and |x^ are related to the frequency of the exploring radio -wave, the elec- 
tfon density in the plasma and the strength and direction of the earth' s field 
through the classical Appleton -Hart ree magneto -ionic theory. Throughout the 
work described here, the calculation of the group refractive indexes (and of 
G 0 and G x ) were made using the Appleton-Hartree equations in the formulation 
given by Becker (ref. 1 6 ). The relevant formulas are summarized in appendix D. 

The electron density for each Ah-j_ is computed from 


Nj_ = 1.2388 x10 4 f|( Ordinary ray) 

iq = 1.2388 x10 4 f^ 2 ( Extraordinary ray). 


(18) 


The use of the single polynomial procedure is limited to cases where the 
N(h) profile is monotonic and the h'(f) curve contains no sharp "cusps." 


The Overlapping Polynomial Method 


The overlapping polynomial procedure described here is based on the work 
of Titheridge (ref. 15) who has applied this technique in the reduction of 
ground-based ionograms. Whereas the single -polynomial method cannot be used 
for reducing ionograms having cusps and other sharp irregularities, the over- 
lapping polynomial method is quite generally applicable. 


Instead of the entire N(h) profile being represented by a single poly- 
nomial, a number of polynomials are used over a series of small frequency 
ranges. Thus, when the real depth is calculated at any frequency, a polyno- 
mial to represent the exact shape near this frequency only is required, and 
virtual depths outside the particular range are not used. The polynomial is 
chosen to give the correct virtual depths at frequencies just above and below 
as well as at the frequency under consideration. The polynomial is also 


joined to an already determined part of 
the curve, and so a smooth continuation 
of the real depth curve with the correct 
gradient is assured. 

In order to illustrate this, sup- 
pose that real depths have been deter- 
mined at frequencies fi, f 2 > . .«,fi_i> 
and that it is desired to compute the 
real depth corresponding to frequency 
f^ (see fig. 2). A fourth-degree poly- 
nomial 


Ah = y a j%^ U9) 

j=° 



f N (fih) KNOWN I f„ (Ah) TO BE COMPUTED 


Fig. 2.- Schematic diagram illustrating the step 
by step and overlapping procedures by which a 
profile is deduced from Ordinary ray ionogram 
data by the overlapping polynomial technique. 
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will be found which smoothly connects the point at f ± on the fjj(Ah) curve 
to the points at fj__ 1 and fj_„ 2 , and such that the correct values for h 1 -*■ 
are maintained at and also at fj^+i- The constants a 0 , a 1; a 3 and 

a 4 are to be computed, following which Ahq is calculated by substitution 
of fj_ for fjj in equation (19)* The five equations to be solved for the 
a ! s are 


Ah 


L m 


4 


J=o 




for m = i - 2 and i - 1 


> ( 20 ) 


r f m , j -i 

J ^ 0 ( f :ni> f N) f N df N > for m=i-l,i,i + l 


0=1 


f i-2 


The group retardations are given by the integral 

'Ahi _ 2 


pAh m 

Ci- 2 = / do(fm.%) d (^h) 

9 v -' Ah • _ 


or 



-2 




( 21 ) 


The numerical treatment of this integral, as well as those in equations (20), 
will he discussed later. 

Special consideration is needed for the last point (point n), since 
point n + 1 does not exist. A third -degree polynomial 


3 



0=0 


( 22 ) 


is used, and the four equations to he solved for Oq, a x , a 2 and a 3 are 
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The group retardations H* are again defined as in equation (21). 

The required two starting values for the overlapping polynomial procedure 
are determined ty use of a lamination method (see following section and table I) 
suitably modified for the topside application. 

The evaluation of the integrals in equations (20), (21) and (23) is facil- 
itated by the introduction of the transformation fjj = f m sin Q (fjj = f^ sin 6 
in the case of the Extraordinary trace). The integrals in equations (20) and 
(23) become 


. T*/ 2 i-i 

l i f N _ “N " J f m J . °o< f m ^ e ) 48 

J sm 1 (f i _ 2 /f m ) 


The calcination of by use of the function G 0 is discussed in appendix D. 

The integral in equation (21) is handled in a similar manner. First, however, 
it is necessary to evaluate the derivative d(Ah/d%) in the frequency range 
of f Nv to . Since Ah has already been calculated in this range, this 
can be done by means of a simple linear approximation between given frequency 
values. Equation (21) then becomes 


Ah k - AhfrA r k , . 

„ *k - *k-i ) j f„ df N 


for m = i - 1, i, and 141 


,, V r /Ah k - Ah k .A r sln " 

‘ A Lv f k - f k-l / m j. in- 1 


-sin- 1 (f k /f m ) 

sin~ 1 (f k _ 1 /f m ) 


0 o d<9 
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(Note that f 0 = % v and Aho = 0.) Equations (19) through (25) apply for tjie 
Extraordinary trace if fj_, p6 an ^ G 0 are replaced by f p, and G x , 

respectively (see appendix D) . 

In the case of Extraordinary trace data, a correction is made, following 
the calculation of the real depths, to adjust for the fact that a constant 
gyrofrequency (namely, fHv) has been used throughout. This correction con- 
sists in approximating fjj-j_ at each value of Ahj_(km) by the inverse cube 
relation 


c Hi 


= ft 


7370. 0 


Hv V 7370.0 - Ah- 


(26) 


I 

and then recalculating fj_ using the relation 



■^xi%i 


(27) 


The electron density is then computed by the use of the formulas 


Nj_ = 1.2388x10 4 f^ 2 , Ordinary trace 


= 1.2388x10 4 fj_ 2 , Extraordinary trace 


(26) 


As noted earlier, the overlapping polynomial process is started by use 
of the linear lamination method. The program is arranged so that any number 
of points and, in fact, the entire profile, may be calculated by the lamina- 
tion method. In practice, a three -point start is adequate for the overlapping 
polynomial program. It is occasionally of interest, however, to use the lami- 
nation method throughout. 


The Linear Lamination Method 

The linear lamination technique for the calculation of N(h) profiles 
from h T (f) curves (refs. 17 and 18) is based on the inversion of the integral 
equation (9) with cp = f^ so that, for example, for the Ordinary ray 



(29) 
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He:re & Ah/dfjj- is assumed constant over each lamination. An outline of the 
method as applied to the topside case is given here, based on a simple modi- 
fication to the procedure for bottomside ionograms described by Thomas and 
Vickers (ref. 19) • It may be noted that certain notation, namely the indexing 
of quantities, is in agreement with that used by Thomas and Vickers and is 
slightly different from that in the foregoing sections. 


The true and virtual heights are assumed to be related by the following 
set of equations : 


h Q = A 0 o Ah o = 0 

hi - A 10 Aho + An Ahi 

h2 = Ago Aho + A21 Ahi + A22 Ah.2 

) (30) 

h^ = A no Ah 0 + ^ni Ah^ + . . . + A^ Ah^ 




Since Ah Q = h^ = 0, as noted above, the first of equations (30) drops out, 
and terms multiplied by Ah,-, in the subsequent equations are zero. The 
resulting set can be written as the matrix equation 


[h 1 ] = [A] [Ah] 


(31) 


The equations represented by equation (31) solved one at a time, lead to 

[Ah] = [B] [h* ] (32) 


which yields the desired values of Ah(fjj) for the given h’(f). 
equation, the elements of [B] are 


In the above 


1 


®nn = a * 

^nn 


B- 


■nn 

A 


■run 


nm - ~ a y 

^nn 


®nm - ® y 


m = n 


m < n 


m > nj 


(33) 
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The values for Ahi, Ah 2 , . 
tion 


. . can he found In succession by use of the re\a. 


II 

= ^ (®n+i,m ^m) + ®n+i,n+i^-n+i 


It is assumed in the lamination technique that d Ah/dfjj is constant 
between adjacent points, so that 


d Ah _ ~ ^m-i 

df N " f m - f m-i 


= const. 


Equation (29) then becomes 


li 

K = Y ( ^ m ' 


where 


W = J f ^ (% - fm-i) ’ m - n 

■ L m-i 


= 0 , m > n J 

(in the case of Mu, f Q = fjj v .) Equation ( 36 ) may be written as the matrix 
equation 

[h>] = [M][D][Ah] (35 


where 


tg 



Wow equation (38) is equivalent to equation (31) and 

[A] = [M][D] (AO) 

Thus, the elements of [A] can he determined from equation (40), following 
which the elements of [B] can he computed using equations (33)- The real 
depths Ahi, • • • are then readily calculated from equation (34) . 

It was pointed out in the previous section that the overlapping polyno- 
mial program employs the lamination method to obtain starting values. The 
entire computation may he carried out using the lamination technique if spe- 
cial control cards are used in this program (see appendix B) . 

The Contiguous Parabola Method 

A second-order polynomial lamination technique for converting topside 
ionograms to N(h) profiles has been described by Doupnik (ref. 20). In prin- 
ciple , the method is similar to the overlapping polynomial technique except 
for two important differences: the polynomial is of second order (i.e., par- 

abolic), and the laminations do not overlap but are contiguous. It is believed 
that either the overlapping polynomial or contiguous parabola programs can be 
applied to topside profiles with equivalent accuracies for the same number of 
scaled h'(f) points » Each can be used for a wider class of ionograms than 
can the single polynomial or lamination programs and each gives higher accu- 
racy than the lamination technique for a given number of h'(f) data points. 

For mathematical detail concerning the contiguous parabola program the 
reader is referred to the report by Doupnik (ref. 20) who has supplied Ames 
Research Center with a "parabolic in In fjj" digital computer program which 
has been modified to accept input data in the standard Ames format (appen- 
dix B). It should be noted that the series of parabolas open along the real 
height axis and hence do not represent "parabolic layers." Parabolas in fjj- 
and in In fjj can be used and each has certain advantages, as discussed by 
Doupnik (ref. 20). The output format for the contiguous parabola program is 
described in appendix C. Some tests of the accuracy of the method are 
included in the next section. 

DISCUSSION OF RESULTS, SELECTED EXAMPLES 

In this section, a description is given of tests made to determine the 
accuracy of the confutations described in the previous section. 

Examples are presented of a number of h'(f) curves used as input data 
to the digital computer. The corresponding N(h) profiles deduced using the 
different techniques described can then be coopared and the relative accuracy 
of each technique assessed. Other factors affecting the suitability of a par- 
ticular technique when applied to routine ionogram reduction are also reviewed 
briefly. An attempt is made to estimate the over -all accuracy which can be 
attained with scaling and reduction techniques that are sufficiently short to 
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Fig. 3*- The accuracy to which an assumed expo- 
nential topside electron density profile can 
be recovered, from an analysis of the corre- 
sponding hx(f) curve, by the single polyno- 
mial technique. 



Fig. 4. - The accuracy to which an assumed expo- 
nential topside electron density profile can 
be recovered, from an analysis of the corre- 
sponding -hi(f) curve, by the contiguous 
parabola technique. 



Fig. 5*- The accuracy to which an assumed expo- 
nential topside electron density profile can 
be recovered, from an analysis of the corre- 
sponding h x (f) curve, by the overlapping 
polynomial technique. 


enable large numbers of ionograms to be 
reduced to N(h) profiles in an econom- 
ical and rapid manner. 

In order to check the accuracy of 
the computations, two kinds of tests 
were devised. In the first of these 
(described below under the heading No- 
Field Case) a model ionosphere is 
assumed such that a simple analytical 
expression exists for the virtual depth 
as a function of frequency. This 
expression is used for computing the 
h ! (f) curve to any desired accuracy. 

This provides a test ionogram which may 
be read into the digital computer as a 
data curve containing essentially zero 
scaling errors. The N(h) profile 
deduced by the different computer reduc- 
tion techniques can then be compared 
with the exactly known starting N(h) 
model and hence the absolute accuracy 
of the computations can be assessed. 

In the second kind of test, an 
actual topside ionogram is read and a 
table formed giving values of h* and 
corresponding values of f. The N(h) 
profiles deduced using the various 
reduction techniques together with the 
tabular values of h ! and f are then 
compared. The errors in reading h* 
and f can be estimated by comparing 
the N(h) profile deduced from a number 
of independent readings of a given 
ionogram. Further details are given in 
the following paragraphs . 

No -Field Case 

When there is no magnetic field, 
it is possible to choose realistic ana- 
lytic approximations to topside profiles 
which are amenable to exact closed -form 
solutions of equations (10), (ll) and 
(12). One such topside model is an 
exponential decay of plasma density with 
real depth corresponding to the curve 

Ah = 400 In f 

and shown as the upper curve in fig- 
gures 3, 4 and 5. The corresponding 


h'(f) curve is given by 

* 

h' = 400 cosh -1 f 


Corresponding values of h’(f) are pre- 
sented in the first two columns of 
table II and plotted as the lower curve 
in figures 3> 4 and 5* 

Values of h* and f computed from 
the above expression were read into the 
digital computer and the W(h) profile 
deduced using the single polynomial, 
contiguous parabola, and overlapping 
polynomial techniques. The correspond- 
ing computed f(Ah) points for each 
technique are shown as circles in fig- 
ures 3, 4 and 5 from which it can be 
seen that the assumed profile is recov- 
ered to a high degree of accuracy. 

Table II shows the values of Ah com- 
puted with the overlapping polynomial 
technique (column 3) together with the 
exact values. It is seen that the 
f(Z^i) curve is recovered everywhere to 
better than -1 km with a mean differ- 
ence between the computed and exact 
W(h) curves of less than 0.4 hm. 
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FREQUENCY, MC/S 

Fig. 6.- Typical equatorial Alouette topside 
ionogram used as standard test for reduc- 
tion of Alouette N(h) data. The ionogram 
was recorded at Singapore on Alouette pass 
number 69 7 on November 19> 1962, at 08l0 TJT. 
The latitude and longitude of the subsatellite 
point was 11.6° S and 117.2° E, respectively. 
At this time the satellite was at an altitude 
of 1003.2 km, and the electron gyrofrequency 
at the vehicle was 0.8l Mc/s. Strong Ordi- 
nary and Extraordinary traces and a weak 
Z-trace are evident in addition to plasma 
and gyrofrequency resonances. 


Magnetic Field Included 


When the magnetic field strength 
is finite, there is no analytic profile 
model available which is amenable to 
exact closed-form solution as above. 

Hence, it is not possible to check the 
absolute accuracy of a given reduction 
technique or associated computer pro- 
gram. However, the relative accuracies 
may be estimated by comparing the N(h) 
profile deduced by each technique from 
a given set of h' (f) points. For 
this purpose the MSA Topside Working 
Group has chosen a working ionogram 

taken from an actual Alouette I frame (fig. 6). Table III lists digitized 
values of h^(f) and h^(f), as scaled from figure 6 by one representative of 
the Topside Working Group. These values are to be used as standard inputs to 
the various reduction programs. The circles and squares in figure 7 are plots 



Fig. 7.- Plot of h£(f) and hy_(f ) for test iono- 
gram of figure 6. See table III for values. 
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of these tabulated points; smooth lines join them. Each member organization 
of the Topside Working Group computed an W(h) profile using the data of ' 
table III. Different reduction techniques were employed and a smooth com- 
posite profile based on all the resulting N(h) data as drawn by the Topside 
Working Group is shown in figure 8 and tabulated in table IV. 

The h 1 (f) data of table III were used as inputs to the Ames Research 
Center reduction programs. The resulting N(h) profiles, for the linear lam- 
ination, contiguous parabola and overlapping polynomial programs are shown in 


figures 9^ 10 and 11, respectively. In 



Fig. 8.- "Composite" Topside Working Group 
electron density profile resulting from 
reduction of the tabulated values of the 
test ionogram by at least eight different 
computer programs as provided by at least 
four different research groups. The graph 
is a "mean" curve through points having a 
typical scatter of ±1 to 2 km. 



Fig. 10.- Test of NASA Ames Research Center tech- 
niques for the reduction of h*(f) curves to 
N(h) profiles (contiguous parabola technique - 
see table I). The input data were the values 
in table III. An inverse cube height variation 
is allowed for in the E trace analysis; con- 
stant fft is assumed in the 0 trace analysis. 


each of these figures the squares 



Fig. 9.- Test of NASA Ames Research Center tech- 
nique for the reduction of h’(f) curves to 
N(h) profiles (lamination method - see 
table I). The input data were the values in 
table III. An inverse cube height variation 
is allowed for in the E trace analysis; con- 
stant % is assumed in the 0 trace analysis 



Fig. 11.- The accuracy with which an overlapping 
polynomial technique produces results compa- 
rable with the "composite" profile of fig- 
ure 8 (continuous curve). The h’(f) values 
of table III were used as input data. An 
inverse cube height variation is allowed for 
in the E trace analysis; constant f H is 
assumed in the 0 trace analysis. 
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represent the profile deduced from the tabulated Extraordinary trace data and 
the circles represent the profile deduced from the tabulated Ordinary trace 
data. The smooth curve is the Topside Working Group’s composite profile of 
figure 8. In all cases, the deduced N(h) profile agrees well with the 
Topside Working Group's curve, with the contiguous parabola and overlapping 
polynomial techniques giving a more accurate profile than the linear lamina- 
tion method. Numerical values of % and Ah deduced from the Extraordinary 
and Ordinary traces, using the overlapping polynomial program, are presented 
in table V. 


Scaling Accuracies 

The main errors in ionogram-profile 
reduction arise from inaccuracies in 
ionogram scaling. Since the basic prob- 
lem concerns the rapid, but accurate, 
reduction of large quantities of topside 
ionograms electronic digitizing machines 
rather than manual scaling techniques, 
have been primarily employed. The 
squares in figure 12 show the Extraor- 
dinary trace of the ionogram of figure 6 
scaled using a Benson-Lehner OSCAR-F 
electronic digitizer. The smooth curve 
of figure 7, representing tabulations 
of the same trace by the Topside Working 
Group, is superimposed. It is clear 
that errors of ~ 10 km in h' can arise 
even when the curves are carefully, but 
fairly rapidly, scaled. Thus random 
variations of ~ 4 km are likely to occur 
in the N(h) profile as a result of 
scaling errors. Figure 13 compares the 
deduced profile using the digitized 
h* (f ) data of figure 12 as inputs to 
the overlapping polynomial program 
(squares) with the profile deduced with 
the same program using the h’ (f ) data 
of table III as inputs (smooth curve). 



Fig. 12.- Comparison of ionogram scaling produced 
"by two independent groups. 



Fig. 13 •- Two profiles computed from the test 
ionogram scaled by two independent groups 
by an overlapping polynomial program. 


Variation of Magnetic Field Strength With Height 

It is important to note, that, especially for the topside ionosphere, it 
is necessary in the computations to allow for the variation of the magnetic 
field strength with height, particularly for Extraordinary trace analysis. 
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Fig. 14. - Improvement of the computed profile 
when a constant electron gyrofrequency is 
replaced by an inverse -cube variation of 
gyrofrequency with depth below the Alouette 
satellite. In both cases the redaction of 
the E trace was carried out by the over- 
lapping polynomial method. 


The circles in figure 14 represent the 
profile deduced from the Ordinary trace 
h'(f) test ionogram (table III) assum- 
ing a constant gyrofrequency with 
height. The squares represent the pro- 
file deduced from the Extraordinary 
trace (table III) again assuming a 
constant gyrofrequency. The smooth 
curve represents the profile deduced 
from the same Extraordinary trace when 
the gyrofrequency is assumed to 
increase with depth according to an 
inverse -cube law. Significantly better 
agreement between the profiles deduced 
from the Ordinary and Extraordinary 
trace data is achieved when the height - 
varying gyrofrequency is employed. 


Nighttime Ionograms 


The ionograms so far presented 
(figs. 3 and 6) have shapes typifying 
daytime conditions in the topside ion- 
osphere. At night the appearance of 
the ionograms is different, the h'(f) 
curve showing a sharp cusp at low fre- 
quencies (fig. 15). Such rapid changes 
of h' with f can severely tax the 
ability of many computer techniques to 
produce a reasonable profile. Clearly, 
the single polynomial technique cannot 
be applied in such cases and neither 
can a lamination technique using large 
values of lamination depths. The con- 
tiguous parabola and overlapping poly- 
nomial techniques are adequate, however, 
and yield satisfactory results. 

Figure 15 illustrates data from an actual Alouette nighttime ionogram. The 
profile deduced by the overlapping polynomial program is also shown in this 
figure and the associated numerical data are presented in table VI. 



Fig. 15*- The manner in which an observed night- 
time ionogram having a deep cusp in the E 
trace corresponds, when computed by the over- 
lapping polynomial technique, to a smoothly 
varying f^(Ah) profile. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. , Mar. 8, 1965 
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APPENDIX A 


INSTRUCTIONS FOR SCALING TOPSIDE h’(f) RECORDS 
INTRODUCTORY COMMENTS 


The purpose of this appendix is to document detailed specifications for 
the scaling of Alouette I ionograms as data inputs to the N(h) digital com- 
puter programs described in the text. In general, these specifications apply 
for the contiguous parabola and overlapping polynomial programs. Minor addi- 
tional considerations pertinent to the lamination and single polynomial pro- 
grams are discussed at the conclusion of the appendix. 

"Scaling" (or "reading") the curve is concerned with the representation 
of an Extraordinary and/or Ordinary topside ionogram hy a finite number of 
pairs of h’ and f values measured at suitable intervals along the trace. 

The accuracy and number of scaled points required for confuting accurate pro- 
files with the techniques available are considered. In practice, the Extraor- 
dinary trace is usually scaled since it is the more frequenty observed in 
its entirety. Hence, reference will be made to that trace only; scaling pro- 
cedures for the Ordinary trace are similar. 


FORMATS 


As described elsewhere in this report, two distinct computer programs 
are utilized for converting ionograms to profiles . One program, described 
in reference 10 and designated here as the Ny program, requires, as iono- 
gram data input, the frequency, f^, at which the Extraordinary trace has 
zero virtual depth, that is, b^(f ) = 0, together with the time of occurrence 
of f^ measured to ±0.5 second. The input ionogram data format for this 
program is shown in table VH. In addition -to the time and f xv , auxiliary 
data required are the pass number, date and telemetry station. Table VIII 
illustrates the positional data input format for the Nv program. This table 
contains the longitude, latitude and height of the satellite for the given 
pass and date at times 1 minute apart such that the time of f^ given in 
table VII falls between these two times. Whether or not the satellite is in 
sunlight at the time stipulated is indicated by an asterisk or a blank in 
column bl; for convenience the planetary index of magnetic activity, K p , at 
the time is also listed. All positional data, in the format of table VIII, 
is punched by hand from "world maps" defining the Alouette I orbit. Table IX 
illustrates the output format for this N v program. Several, but not all, 
of the output parameters are required as inputs to the N(h) programs. 

Tables X and XI (appendix B) show the identification card input format 
and parameter card input format for the N(h) programs. The parameters com- 
mon to tables IX, X and XI are pass number, date, time, dip angle. 
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gyrofrequency at the vehicle, f Hv and f xv . It should be noted that the 
quantities dip angle and f^y may be entered in the N(h) input format only 
after the N v program has been run. Table XII illustrates the frequency- 
virtual depth data input format which contains all h' and f pairs for a 
given ionogram excluding the f xv value. The reader is referred to refer- 
ence 10 and appendix B herein for further necessary detail on input card 
formats. 


IONOGRAM QUALITY AND VARIABILITY 

Figures 1 6 and 17 illustrate "good" Alouette ionograms inasmuch as the 
Extraordinary trace has no large gaps, is narrow in width (< 0.1 Mc/s in 
thickness at small h 1 and < 20 virtual km in thickness at greater h 1 ), and 
has fairly sharp contrast. These ionograms are rather easily scaled. Fig- 
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Fig. 16.- "Good" Alouette I ionogram illustrating 
typical daytime traces and the satellite 
trace. 
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Fig. 17.- "Good" Alouette I ionogram illustrating 
typical nighttime Extraordinary and satellite 
traces . 


ure 16 , typical of daytime quiet con- 
ditions, shows h^. to be a smooth, 
monotonically increasing, single -valued 
function of fx- On the other hand, 
figure 17 is typical of nighttime quiet 
conditions where f x is a multivalued 
function of iu. due to the evident 
cusp (large group retardation) at fre- 
quencies only slightly higher than 
^xv Considerable care must be exer- 
cised in the scaling of the latter ion- 
ogram and additional h 1 and f pairs 
may be required for equivalent accuracy 
(see below). Figure 18 illustrates a 
"poor" ionogram typical of "spread" 
ionospheric conditions. Clearly, any 
attempt to scale the latter would lead 



f, MC/S 


Fig. l8.- "Poor," hut scalable, Alouette I iono- 
gram illustrating typical disturbed ionospheric 
conditions . 
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only to a rough, first-order estimate of the real height. The ionograms which 
. will he encountered in analyses of Alouette I records will range in quality 
and variability between those illustrated in figures 1 6 to l8. It should be 
noted from these figures that the satellite trace (equivalent to the ground 
pulse in ground-based sounders) provides the zero virtual depth measurement 
and must therefore also be of "good" quality. 


GENERAL SCALING PROCEDURE 


The basic requirement is to tabulate h^ at a set of discrete frequen- 
cies for the Extraordinary trace from f^ to f x F 2 when the trace extends 
that far, or to the highest observable f x . The Extraordinary trace is always 
the one that begins and ends farthest to the right on a given ionogram. 

All h* and f points should be scaled at the leading edge of the trace 
(see figs. 1 6 to 18 ). 

The zero -depth reference is always the leading, that is, top, edge of 
the satellite trace. The zero height marker should coincide with the leading 
edge of the satellite trace. Frequency calibration is afforded by the fre- 
quency markers. Linear interpolation between height markers and frequency 
markers can be used. 

In general, scaling should only be attempted, and will only be reliable, 
when the leading edge is fairly well defined and when no large gaps exist in 
the trace (see figs. 16 and 17). Ionograms with gaps of less than 0.5 Mc/s, 
or even in some cases up to 1 Mc/s, may be scaled by simply reading a smoothly 
drawn curve through the gap bearing in mind the shape of the trace on similar 
days or on other ionograms within the pass. 

When the h' (f ) trace shows signs of severe scattering or masking by 
resonances, spurious responses, spread F, overlapping of Ordinary and Extraor - 
dinary traces, etc., it may not be possible to scale a given ionogram. If the 
masking is localized, particularly near f x F 2 , it may still be worthwhile to 
scale the ionogram by drawing a smooth curve across the locally masked area 
bearing in mind the shape of the trace on similar days or on other ionograms 
within the pass. 

From an ionospheric physics viewpoint, very interesting effects are often 
prevalent on ionograms recorded during disturbed, rather than quiet, condi- 
tions. Hence, with some sacrifice in reliability, it may be advantageous to 
scale certain types of spread ionograms. An example is given in figure l8. 

It should be noted that this ionogram has an easily identified Extraordinary 
trace but has considerably more spread than the corresponding trace of fig- 
ures 16 and 17. The trace has a rather well-defined trailing edge from 0 
down to about 800 virtual km and a comparatively well-defined leading edge 
thereafter. One can therefore construct a scalable trace by drawing a smooth 
curve, from 0 to 800 virtual km, vertically above the trailing edge by an 
amount roughly equivalent to the width of the satellite trace, approximately 
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20 virtual km in such a manner that it smoothly joins the leading edge at 
greater depths. Such data should he marked with a sign to indicate this spe-- 
cial procedure has had to he used. 

SPECIFIC SCALING DETAILS 


There is no absolute criterion for the number and density of scaled 
points. It is necessary that a minimum number of scaled points, meeting the 
criteria below, be used. Overscaling decreases the number of ionograms which 
can be read, per unit time, with no increase in accuracy. A minimum of 8 to 
10 points is usually needed, rarely more than 15 to 20. 

For traces of the type shown in figure l6 the leading edge is scaled at 
about 0.05 to 0.1 Mc/s intervals for the first few points. The remainder of 
the curve is scaled at intervals of about 0.5 to 1 Mc/s. For an accurate 



Fig. 19.- typical daytime and nighttime sketched 
Extraordinary traces with suggested scaling 
locations . 

about 10 pairs of h* and f values ar 
"low" cusps with as many as 20 require 


N(h) profile, these remaining values 
should be very carefully scaled. A 
typical example is given in curve a of 
figure 19. A total of not more than 
10 pairs of values of h' and f should 
be sufficient. 

For traces of the type shown in 
figure 17 it is necessary to scale in 
steps of between 0.01 and 0.05 Mc/s 
from f xv down to and around the cusp, 
so that at least f or 8 points, includ- 
ing that at the cusp, are recorded in 
this interval. Following this, scaling 
is done in intervals of 0.2 to 1 Mc/s. 
Examples are given in curve b of fig- 
ure 19 and in figure 15 . A minimum of 
; normally required for ionograms having 
1 for more extreme cases. 


It is always necessary to ensure that successive frequency readings 
increase monotonically (assuming that the trace is scaled in the direction 
from f xv to f x F 2 ) . 

Irregularities with a scale of < 0.1 Mc/s or < 5 virtual km should be 
smoothed out unless otherwise advised (e.g., for special studies). 

The input formats require h^. to be rounded off to the nearest km and 
f x to the nearest 0.01 Mc/s. 


LAMINATION PROGRAM 


In principle, it is possible to achieve equivalent accuracies with the 
lamination technique as compared with the contiguous parabola or overlapping 
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polynomial techniques. However, this requires a significant increase in the 
pumber of pairs of h' and f values (generally two to three times as many, 
Thomas and Vickers (ref. 19)) otherwise, the above scaling procedures apply. 


SINGLE POLYNOMIAL' PROGRAM 


Because of the relatively large number of pairs of h’ and f values 
required in the lamination .program, the single polynomial and overlapping 
polynomial techniques (ref. 15) and contiguous parabola technique (Doupnik, 
ref. 20) were devised. Although the latter three techniques thus represent 
a significant improvement over the classical lamination approach, the single 
polynomial program suffers from the fact that a maximum of eight to nine pairs 
of h' and f values are useable (Thomas and Westover, ref. 14). Hence, this 
program is applicable to ionograms only of the class illustrated in figures l6 
and l8. Ionograms showing cusps (as fig. 17) should never be reduced using 
the single polynomial program. Otherwise, the above scaling procedures apply. 
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APPENDIX B 


INPUT FORMATS 


The observed topside ionogram data (the h'(f) curves) are read into the 
digital computer on punched cards , the same format being used for each of the 
different programs. Special control cards are, however, required for select- 
ing certain program options. 

The standard input consists of three types of cards. These are an iden- 
tification card, a card containing parameters, and a group of cards containing 
frequency and virtual height pairs, four pairs per card. Formats for these 
three types of input cards are shown in tables X, XI and XII. 

It should be noted (table XI) that if the decimal point is included, the 
quantities I Y , % Y and f Y (or f XY ) may be placed anywhere within the speci- 
fied fields. The decimal point may be deleted, in which case it is understood 
to lie between the second and third positions from the right in the I Y field, 
and between the third and fourth positions from the right in the fjjy and f Y 
(or f XY ) fields. 

The frequency and virtual depth data (table XII) may be punched with any 
format whatever, in the specified fields, if the decimal point is included. 

If the decimal point is deleted, it is understood to lie between the second 
and third positions from the right in the frequency fields, and to the right 
of the right-most position in the virtual depth fields. It should be pointed 
out that the first frequency and virtual depth pair on the first of a set of 
data cards is not the reading corresponding to zero virtual depth, but rather 
is the first reading at a finite virtual depth. 

The single polynomial program and the contiguous parabola program will 
accept standard data as described above. A summary table of interpolated 
N(Ah) values, which is routinely printed out by the single polynomial program, 
may be deleted if a "1" is punched in column . 72 of the ionogram identification 
card. (This summary table is never prepared by the contiguous parabola pro- 
gram. ) 

The overlapping polynomial program is normally used with a three -point 
linear lamination method start, and interpolated tables of N versus Ah and 
scale height versus Ah are routinely saved. This procedure uses the stand- 
ard data, as described in the foregoing paragraphs. If it is desired to use 
other than a three -point lamination start (or to compute the entire ionogram 
with the lamination method), or if it is desired to delete the two interpo- 
lated tables, then a modified version of the overlapping polynomial program 
is used, and a control card of the format described in table XIII is placed 
behind the data for each ionogram to be computed. 

In all cases , a blank card is placed behind the last card in a run. 
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APPENDIX C 


OUTPUT FORMATS 


Detailed printouts are made for each ionogram. Examples are shown for 
an Extraordinary trace computation using the overlapping polynomial program 
and the contiguous parabola program in tables XIV and XV, respectively. It 
should be noted that pass number, date, time, dip angle Iv> gyrofrequency at 
the vehicle fg v and Extraordinary ray frequency at the vehicle f xv , appear 
in the heading. Then columns of frequency, real and virtual depth and elec- 
tron density are printed. The column headings for the Extraordinary ray 
n detailed output” are defined below. 


Quantity 


Overlapping polynomial 


Contiguous parabola 


h 1 H PRIME H 

f x EX F 

f 1 F PRIME FN 

Ah DELTA H DELTA H 

N ELECTRON DENSITY DENSITY 

Real height - - - Z 

Gyrofrequency - - - FH 


The irwo columns labelled A and B in table XV are intermediate quantities and 
have no physical meaning. 

In the case of an Ordinary trace computation, the column for f x (ex) 
does not appear, and f' (F PRIME) is replaced by f (F), in the overlapping 
polynomial program detailed output. Column headings are not changed in the 
contiguous parabola program; however, and the quantity f apears in both 
columns F and FN of the detailed printout. 


A s umm ary interpolated table of electron density versus real depth is 
available as output from the single polynomial and overlapping polynomial 
programs. An example is shown in table XVI. In addition, a summary interpo- 
lated table of scale height versus real depth is available as output from the 
overlapping polynomial program. The scale height is calculated by means of 
the formula 


2AhiNi 

Scale height = t: 77 

i'i+i - «i-i 


where here, Ahq = Ah(%) - Ah(Ni_ x ) = Ah(N i+1 ) - Ah(Nq) = 50 km. An example 
of this output is shown in table XVII. 
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Approximate computing times have been determined on an IBM 709^ computer 
as follows : 

(a) Single polynomial with N(Ah) summary table, a 7-reading ionogram, 
0.015 minute. 

(b) Overlapping polynomial with N(Ah) summary table, a 14 -point iono- 
gram, 0.040 minute. 

(c) Contiguous parabola program, a 14 -point ionogram, 0.080 minute. 
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APPENDIX D 


CALCULATION OF GROUP REFRACTIVE INDICES 


The Appleton-Hartree magneto -ionic equations in the formulation given by- 
Becker (ref. 16) are used here to compute the group refractive indices 
and jx^. The following definitions are made; first. 


Ordinary ray 


Extraor dinary ray 


Uq = Go/t 0 

P x - G x /t x 

t o = 1 - 

4 - 1 - (r; 

y Q = f H / f i 

y x = %/ f j 

X o = 1 - *o 

^ - 1 - 1; 


Then, for the Ordinary ray, 

Xo tan 2 I 


to 

G Q = — 
° n o 


1 + 


W 


1 + X r 


Jl + 7 ^ 1 + Jl + 7t 


(Dl) 


where 


7 = 


jfc tan 2 I 
y^ cos 2 ] 


M = 1 + t‘ 


2 2 tan 2 I 


!+>/! + 7tc 


1 + 


no\ 2 

(W 


2 tan 2 I 

1 + Jl + /to 




1 + t' 


2 tan 2 I 

° 1 + Jl + TtJ 


(D2) 
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and for the Extraordinary ray. 


p _ 

x " 


1 + 


N 2 


1 + l 


2| (1 + j)(l + X)> 


. 1 + \/ 1 + «7 1 + cl | 


(D3) 


where 


a = 


4 sin 2 I 
(l + sin 2 I) 2 


(3 = 


2 sin 2 I 
1 + sin 2 I 


\ = 2A + 


A = t 2 , 

X y ' 


= I y x C / 4 + y x - y^) 


X = (l - y)x. 


X 


■h 




N=l-y+A- 


J&L 


1 + \/l + 


2 . 1 + sll + a* 

Hx - y^ n 


(D4) 


Some limiting conditions are: 


lim 

t 0 -*o 


n r 




1 

cos I 


lim G q = 1 
t o - " 1 



•N 

2 

sin 2 I + 1 


lim G x = 1 

V* 1 


(D5) 


(d6) 
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TABLE I.- TABLE OF REFERENCES. [THE NASA AMES RESEARCH CENTER DIGITAL 
COMPUTER PROGRAMS FOR h’(f) - N(h) REDUCTION ARE DIRECTLY BASED ON 
THE TECHNIQUES DESCRIBED IN THE MIDDLE COLUMN] 


Method 


Direct references 


Important related 
references 


Single polynomial Titheridge , 1961 (ref. 15) Unz, 1961 (ref. 21) 


Overlapping 

polynomial 



Contiguous 

parabola 


General reviews 


Thomas, Long and Westover, 
1963 (ref. 13) 

Thomas and Westover, 1963 
(ref. 14) 


Titheridge, 1961 (ref. 15) 


Budden, 1954 (ref. 17) 
King, 1954 (ref. 18) 


Thomas and Vickers, 1959 
(ref. 19) 


Doupnik, 1963 (ref. 20) 


Knecht, VanZandt and 
Watts, 1962 (ref. 22) 




Paul and Wright, 1963 
(ref. 23) 

Fit zenre iter and Blumle, 
1964 (ref. 24) 

Jackson, 1956, 1963* 

( refs . 25 , 26 ) 


Doupnik and Schmerling, 
1965 (ref. 27) 


Thomas, 1959 (ref. 28) 


Appleton -Hart ree 
magneto -ionic 
theory 


Becker, i960 (ref. l6) 


Ratcliffe, 1959 (ref. 29)1 






















_ TABLE II.- COMPARISON OF EXPONENTIAL MODEL PROFILE AND THE PROFILE DEDUCED 
FROM THE CORRESPONDING h’(f) CURVE USING THE OVERLAPPING POLYNOMIAL 
METHOD 


h', 

km 

Mc/s 

Ah 

km 

y 

Confuted. 

Exact 

0.0 

1.00 

0.0 

0.0 

177 A 

1.10 

37-5 

38.12 

21+8.9 

1.20 

72.0 

72.93 

302.6 

1.30 

103-9 

104.94 

1 346.8 

1.40 

134.0 

134.59 

385.0 

1.50 

161.6 

162.19 

463 '. 5 

1-75 

223.5 

223.85 

526.8 

2.00 

276.9 

277.26 

626.7 

2.50 

366.0 

366.52 

705.1 

3-00 

439.1 

439.44 

769.9 

3-50 

500.6 

501.10 

825.4 

4.00 

553-9 

554.52 

878.9 

4.50 

602.1 

601.63 

917.0 

5.00 

644.2 

643.78 

955-8 

5.50 

681.8 

681.90 

991.2 

6.00 

716.6 

716.70 

1024.0 

6.50 

748.6 

748.72 

1054.0 

7.00 

778.3 

778.36 

1081.0 

7.50 

805.7 

805.96 

1107.0 

8.00 

831.4 

831.78 
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TABLE III.- TOPSIDE WORKING GROUP TEST IONOGRAM 


Telemetry station: Singapore Alouette altitude: 1003.2 km 

Pass: 697 (19 Nov. 1962) Subsatellite longitude: 117.2° E 

Frame: 28(0810 UT) Subsatellite latitude: 11.6° S 

fjj = 0.8l Mc/s; Dip angle = 4l° 

(h^ and h^. are Ordinary and Extraordinary ray virtual depths (km), and 
f and fy . are Ordinary and Extraordinary ray frequencies of reflection 
(Mc/s)) 


Ordinary trace 

Extraordinary 

trace 


f 


f 

4 

f X 

hi 

f x 

0 

1.63 

730 

7*50 

0 

2.08 

700 

5.50 

230 

1.70 

750 

8.00 

175 

2.10 

705 

6.00 

355 

1.80 

765 

8.50 

340 

2.20 

710 

6.50 

415 

1.90 

795 

9.00 

430 

2.30 

720 

7.00 

46o 

2.00 

805 

9.10 

480 

2.40 

735 

7.50 

490 

2.10 

810 

9.20 

525 

2.50 

750 

8.00 

510 

2.20 

820 

9.30 

555 

2.60 

770 

8.50 

530 

2.30 

8 30 

9.40 

580 

2.70 

790 

9.00 

545 

2.40 

840 

9.50 

600 

2.80 

795 

9.10 

560 

2.50 

850 

9.60 

615 

2.90 

800 

9.20 

575 

2.60 

860 

9.70 

630 

3-00 

805 

9.30 

585 

2.70 

870 

9.80 

640 

3.10 

810 

9.40 

595 

2.80 

885 

9.90 

645 

3.20 

815 

9.50 

605 

2.90 

900 ! 

10.00 

655 

3-30 

820 

9.60 

610 

3.00 

920 I 

10.10 

660 

3.40 : 

830 

9.70 

635 

3.50 1 

945 

10.20 

665 

3.50 

835 

9.80 

650 

4.00 

975 

10.30 

670 

3.60 

845 

9.90 

660 

4.50 



670 

3.70 

850 

10.00 

675 

5.00 



675 

3.80 

860 

10.10 

685 

5.50 



680 

3.90 

870 

10.20 

695 

6.00 



680 

4.00 

880 

10.30 

705 

6.50 



690 

4.50 

890 

10.40 

715 

7.00 



695 

5.00 

905 

10.50 






TABLE V. - REDUCTION OF THE TOPSIDE WORKING GROUP TEST IONOGRAM USING THE 

METHOD OF OVERLAPPING POLYNOMIALS 


Ordinary trace 

Extraordinary trace 

h’. 


Ah, 


fx, 

f N> 

Ah, 

km 

Mc/s 

km 

km 

Mc/s 

Mc/s 

km 

0 

1.63 

0 

0 

2.08 

I.63 

0 

230.0 

1.70 

25.7 

175.0 

2.10 

1.64 

9.0 

355.0 

1.80 

62.8 

340.0 

2.20 

1.74 

45.1 

415.0 

1.90 

96.7 

430.0 

2.30 

1.84 

78.9 

460.0 

2.00 

127.5 

480.0 

2.40 

1-93 

111.0 

510.0 

2.20 

180.6 

525.0 

2.50 

2.03 

138.3 

560.0 

2.50 

244.4 

630.0 

3.00 

2.51 

254.5 

610.0 

3-00 

323.9 

665.0 

3.50 

3.00 

329.8 

635.0 

3.50 

380.3 

680.0 

4.00 

3.49 

383.9 

650.0 

4.00 

421.7 

690.9 

4.50 

3-99 

424.6 

660.0 

4.50 

453-4 

695.0 

5-00 

4.48 

456.4 

675.0 

5-00 

479.8 

700.0 

5.50 

4.98 

481.9 

685.0 

5-50 

502.3 

705.0 

6.00 

5-48 

503.3 

695.O 

6.00 

521.4 

710.0 

6.50 

5.98 

521.5 

705.O 

6.50 

538.4 j 

720.0 

7.00 

6.47 

538.0 

715.0 

7-00 

553-7 | 

735-0 

7-50 

6.97 

554.0 

730.0 

7.50 

568.3 

750.0 

8.00 

7.47 

569.5 

750.0 

8.00 

583.2 

770.0 

8.50 

7-97 

585.0 

765.O 

8.50 

597-5 

790.0 

9.00 

8.46 

600.5 

795-0 

9.00 

612.8 

800.0 

9.20 

8.66 

606.8 

840.0 

9-50 

631.0 

810.0 , 

9.40 

8.86 

613.2 

900.0 

10.00 

653.0 

820.0 

9.60 

9.06 

619.7 


835.0 

9.80 

9.26 

626.5 


1 


850.0 

10.00 

9-46 

633.8 

















TABLE VI.- REDUCTION OF A NIGHTTIME IONOGRAM USING THE METHOD OF 

OVERLAPPING POLYNOMIALS 


Telemetry Station: Stanford- 

Pass: 3901 (12 July 1963) 

Time: 2256 UT 

fjj v = 0.8l Mc/s; Dip angle = 53-78° 


h', 

km 

f x> 

Mc/s 

f N> 

Mc/s 

Ah, 

km 

0 

1.33 

0.83 

0 

468.0 

1.36 

.86 

36.0 

725-0 

1.40 

.89 

83.7 

971-0 

1.43 

.91 

128.0 

IO99.O 

1.48 

.95 

193.8 

1007.0 

1.64 

1.09 

314.3 

876.0 

1.89 

1.33 

401.7 

806.0 

2.22 

1.66 

467.3 

790.0 

2.64 

2.08 

522.0 

803.O 

3.13 

2.57 

569.5 

847-0 

3.67 

3.10 

6l6.2 

935-0 

4.30 

3.75 

673-4 



TABLE VII.- IONOGRAM DATA INPUT FORMAT FOR N v PROGRAM 
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TABLE X. - N(h) IDENTIFICATION CARD FORMAT 


Column 


4-8 


9-19 

20-21 

22-23 

24-25 


Item 

Pass number, integer, right adjusted 
Date, any format whatever (e.g., 10 JULY 63, 


Hours 

Minutes 

Seconds 


Time 


or 7/10/63, etc.) 


Column 


3 


4-6 


7-12 

13-18 

19-24 


TABLE XI.- N(h) PARAMETER CARD FORMAT 


Item 


Control digit 


0 Ordinary ray data 


J 1 Extraordinary ray data 

The number of points read from the ionogram, integer, right 
adjusted 

I v , dip angle at the position of the satellite, degrees (±XX.XX) 

f Hv , gyrofrequency at the position of the satellite, Mc/s (X.XXXX) 

f*v (or f xv ) , the frequency of the Ordinary ray (or Extraordinary 
ray) at the position of the satellite, Mc/s (X-XXXX) 


42 


TABLE XII.- FREQUENCY -VIRTUAL DEPTH DATA CARDS 


Columns 

Item 




2-6 

f (or f x ), Mc/s 

1 reading 

i 


7-12 

h 1 , km 




13-17 

f (or f x ), Mc/s 

1 reading 

i + 

1 

18-23 

h*, km 




24-28 

f (or f x ), Mc/s 

1 reading 

i + 

2 

29-34 

h* , km j 




35-39 

f (or f x ), Mc/s 

» reading 

i + 

3 

40-45 

h ' , km 



1 


TABLE XIII.- SPECIAL CONTROL CARD FORMAT 8, 


Columns Item 

1-5 The number of points to "be computed using the linear lamination 

method, right adjusted. 

6-10 Logical tape upon 'which the two summary tables are to "be written, 

right adjusted. If this field is blank, or zero, then the table 
preparation is deleted. 


a Table XIII is not normally required in routine h T (f)-N(h) reduction, the 
overlapping polynomial method being used with a ” three linear lamination” 
starting procedure. 


l 
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TABLE XIV. - TYPICAL DETAIL OUTPUT PAGE, OVERLAPPING POLYNOMIAL PROGRAM 
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